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I.  INTRODUCTION 


In  a  vide  variety  of  Air  Force  applications,  highly  detailed  information 
about  atoms,  molecules,  and  their  interactions  is  required  [1-3].  This 
information  is  necessary  in  problems  ranging  from  chemical  laser  development  >  to 
the  detection  and  identification  of  rocket  plumes,  to  clustering  and  aerosol 
formations,  non-linear  optics,  electron  beam  technology,  rocket  fuels,  and  even 
to  nuclear  weapons  effects  [1-3]. * 

The  crucial  component  needed  to  understand  molecular  reactions  is  the 
potential  energy  surfaces  that  serve  to  describe  the  attractions  among  the  atoms 
and  molecules  [1].  However,  such  information  is  not  easy  to  obtain.  A  certain 
amount  of  information  about  the  molecular  forces  near  equilibrium  in  a  bound 
molecule  is  available  from  spectroscopy.  Some  information  about  the  potential 
energy  surface  even  in  the  absence  of  binding  can  be  provided  from  crossed 
molecular-beam  experiments.  But,  in  general,  potential  energy  surfaces  are  not 
amendable  to  experimental  determination.  Instead,  other  types  of  observations 
such  as  kinetics  experiments,  coupled  with  very  simple  theoretical  models  for  a 
surface,  are  used  to  infer  pieces  of  information  about  the  parameters  of  the 
model  such  as  what  the  activation  barrier  might  be.  In  addition  there  is  an 
unrelenting  need  for  essential  information  about  other  molecular  properties  like 
excitation  energies,  oscillator  strengths,  moments,  polarizabilities,  etc. 

In  many  cases,  the  most  direct  approach  to  obtaining  accurate  potential 
surfaces  for  molecules,  and  detailed  information  about  their  excited  states, 
vibrational  spectra,  and  a  wealth  of  other  quantities,  is  high  level  ab  initio 
quantum  mechanical  calculations. 

^References  are  on  page  14. 
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However,  more  so  than  in  aost  other  areas,  the  ability  to  provide  reliable 
quantum  mechanical  results  for  increasingly  large  molecules  depends  critically 
on  improved  method  development.  Whereas  supercomputers  can  enable  us  to  make 
much  larger  computations  vith  the  same  old  methods,  the  simultaneous  development 
of  new  methods  can  increase  computational  capability  by  a  further  two  or  three 
orders  of  magnitude.  In  this  regard,  many-body  perturbation  theory  (MBPT) 

[4-15]  and  its  infinite-order  extensions  termed  coupled-cluster  (CC)  methods 
111,16-22]  offer  a  number  of  attractive  features  that  the  more  traditional 
configuration  interaction  approaches  lack  [21]. 

Under  AFOSR  support,  ve  have  established  these  CC/MBPT  theories  as  being 
among  the  aost  accurate  available,  and  have  developed  very  efficient  and 
generally  applicable  computer  codes  to  perform  CC/MBPT  calculations.  Also,  ve 
have  employed  these  methods  for  the  first  time  in  large-scale  ab  initio 
calculations  of  potential  energy  surfaces  [11,21].  The  successes  of  our 
original  work  in  this  effort  have  been  substantial  (see  previous  AFOSR  reports.) 
Tvo  of  our  papers  [11,22]  supported  by  AFOSR  have  been  identified  in  a  book  by 
Schaefer,  Quantum  Chemistry  (Oxford,  1984),  as  being  among  the  149  most 
influential  papers  in  the  50-year  history  of  computational  quantum  chemistry. 
Recently  four  of  our  papers  supported  by  AFOSR  have  been  identified  by  Current 
Contents  as  qualifying  as  science  citation  classics  [7, 10a, 11, 21] ,  and  one  vas 
chosen  to  be  featured  [23].  In  each  case  the  paper  has  been  cited  more  than  tvo 
hundred  times  and  found  to  be  among  the  ten  most  cited  papers  in  the  respective 
journals. 

Following  a  statement  of  research  objectives,  ve  reviev  the  recent 
scientific  advances  ve  have  made  under  AFOSR  support. 


II.  RESEARCH  OBJECTIVES 

The  overall  objective  of  our  continuing  research  program  for  AFOSR  has 
included  the  following: 

(1)  The  development  of  nev,  more  accurate  and  more  efficient  ab  Initio  quantum 
mechanical  methods  based  upon  coupled-duster  (CC)  and  many-body 
perturbation  theory  (MBPT)  for  determining  molecular  properties  and 
potential  energy  surfaces  for  molecular  interactions. 

(2)  The  implementation  of  these  methods  into  highly  efficient,  transportable 
computer  codes,  to  enable  computations  on  molecules  to  be  made  on  an  almost 
routine  basis,  for  a  vide  variety  of  different  properties. 

(3)  The  application  of  these  techniques  to  a  variety  of  problems  that  are  of 
Interest  to  AFOSR,  and  that  serve  to  establish  the  range  of  accuracy  for 
CC/MBPT  methods. 

The  underlying  goal  of  our  effort  is  highly  accurate  first  principle  quantum 
mechanical  predictions  of  molecular  properties.  Such  tools  are  already  making 
possible  the  determination  of  dissociation  energies,  moments,  polarizabilities, 
excitation  energies,  and  vibrational  and  rotational  spectra  for  molecules  to  an 
accuracy  that  is  predictive  and  comparable  to  experiment.  Furthermore,  for 
transient  molecules  as  occur  in  rocket  plumes,  combustion,  flame  chemistry,  and 
interstellar  space  that  are  generally  difficult  to  study  experimentally, 
theoretical  calculations  will  provide  otherwise  inaccessible  information.  In 
addition,  theoretical  work  provides  the  underlying  framevork  to  understand 
innumerable  important  molecular  phenomena  and  to  suggest  previously 
unanticipated  solutions  to  defense  related  problems. 
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III.  NOTABLE  ACCOMPLISHMENTS  (November  1,  1987  -  October  31,  1989) 

A.  For  the  first  time',  contributions  from  "connected"  quadruple  excitations, 

T^,  have  been  included  into  coupled-cluster  theory.  This  allovs  us  to 
define  the  CCSDTQ-1  method  [24].  Comparisons  vith  reference  full  Cl  results 
demonstrate  an  average  error  of  only  0.35  mh  (0.22  kcal/mol).  This  should 
be  compared  vith  the  best  prior  method,  CCSDT,  vhich  ve  presented  a  few 
years  ago  [25],  vhose  error  for  the  same  set  of  examples  is  0.93  ah. 

o 

B.  The  direct  inclusion  of  T^  into  CC  theory  involves  an  n  step,  vhere  n  is 
the  dimension  of  the  basis  set.  This  is  one  pover  of  n  vorse  than  CCSDT. 
Deriving  from  our  recently  proposed  XCC  and  UCC  methods  [26-30],  ve  vere 
able  to  shov  that  the  initial  contribution  of  T^  could  be  introduced  vith  no 
vorse  than  an  n7  step.  Exploiting  this,  ve  proposed  a  highly  efficient  non¬ 
iterative  Inclusion  of  T^  to  be  added  to  any  CCSDT-n  model.  This  Q(CCSDT) 
method  [24]  vas  found  to  very  accurately  reproduce  full  Cl  results  to  a  mean 
error  of  only  0.23  mh.  As  this  is  the  most  efficient  method  that  includes 

7 

T^  (no  step  is  vorse  than  n  )  it  should  find  vide  applicability.  It  also 
offers  significant  improvement  for  the  extreme  test  case  of  Be2  [24]  (see 
Figure  2). 

C.  The  first  report  of  results  for  the  nev  fourth-  and  fifth-order  expectation 
value  XCC  methods  vere  made  [28,29].  Unlike  the  usual  CC  method,  for 

T?  t 

the  energy  XCC  uses  the  ansatz  that  AE><o|(e  He  )g|o>  and  approximations 

^re  obtained  by  energy  variation  to  some  order  in  perturbation  theory. 

XCC(4)  [29]  vas  shovn  to  provide  results  falling  betveen  MBPT(4),  the 
standard  used  in  most  chemical  applications  today,  and  our  CCSDT- 1  method 
[31].  The  fifth-order  version  XCC(5)  along  vith  CCSDTQ-1  and  CCSDT+Q ( CCSDT ) 
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are  the  first  CC  Methods  correct  through  fifth-order.  Barring  convergence 
difficulties,  this  method  vas  found  to  be  quite  accurate  but  not  as  good  as 
the  other  fifth-order  methods.  A  modified  version  called  XCC(5)M  overcame 
the  convergence  problems,  and  had  a  mean  error  of  0.84  ah  compared  to  full 
Cl  [28]. 

D.  The  other  nev  method  ve  investigated  vas  the  unitary,  UCC  method  127].  Using 

t  J  T  T 

t«T-T',  AE«<o|(E  He  >c|o>-<o|(He  ){,|o>.  Energy  variation  leads  to  the 
UCC(n)  equations.  Unlike  other  CC  methods,  UCC  is  guaranteed  to  satisfy  the 
generalized  Hellmann-Feynman  theorem.  This  makes  it  far  easier  to  evaluate 
properties  other  than  the  energy.  Ve  Implemented  UCC(4)  for  analytical 
gradients  vhich  facilitates  the  search  for  minima  and  transition  states  on 
energy  surfaces  and  for  the  determination  force  constants  vhich  leads  to  the 
prediction  of  vibrational  spectra.  UCC(4)  vas  found  to  give  very  good 
structures  and  vibrational  frequencies  for  several  examples  [30]. 

E.  A  byproduct  of  the  above  vas  that  ve  vere  able  to  report  general  results  for 
fifth-order  MBPT  for  the  first  time  [32].  MBPT(5)  constitutes  the  reference 
point  for  all  CC  methods  correct  through  fifth-order.  Ve  reported  results 
for  a  number  of  examples  including  0^  vhich  is  knovn  to  require  a  multi- 
configurational  zeroth-order  description.  Consequently,  a  single  reference 
method  should  shov  large  errors  in  higher-order.  Ve  find  that  E^  amounts  to 
20  mh  and  that  contributions  from  connected  quadruples  is  a  very  large  -16 
mh. 


F.  "One  of  our  other  primary  accomplishments  this  year  has  been  the  continued 
development  of  our  "relaxed  density"  correlated  method  [33-36].  This  nev 
density  is  fundamental  to  many  studies  of  chemical  bonding.  In  our  study  of 
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the  bond-breaking  reaction  CH^  -»  CH  +  H,  ve  used  color  graphics  to  plot  r  I  1, 
correlated  relaxed  densities  to  facilitate  an  understanding  of  correlation 
effects  on  this  bond  breaking  reaction  [37].  'ru 

s 

G.  The  relaxed  density  permits  the  facile  evaluation  of  any  one-electron  e 

property.  This  vas  illustrated  by  our  study  of  numerous  properties  like  /v 

field  gradients,  dipole  and  quadrupole  moments,  and  several  other  properties  e  p 
[35,36]. 

o 

H.  The  relaxed  density  is  also  the  cornerstone  of  all  our  recently  developed 

analytical  gradient  methods.  In  the  last  year,  ve  have  generalized  such  / 

techniques  to  include  the  full  MBPT(4)  method  and  CCD  [36].  [ 

OF 

I.  In  another  major  advance,  ve  presented  the  theory  for  full  CCSDT  analytical  ">1 

gradients  [38]  and  CCSDT-1  analytical  second  derivatives  [39].  The  general 
theory  has  other  methods  such  as  Pople's  so-called  "quadratic  Cl"  as  special 
cases .  z 

ca 

J.  In  a  study  demonstrating  the  beneficial  mutual  interaction  of  theory  vith  ; 

experiment,  an  experiment  of  Zare,  et  al.  [40]  for  the  barrier  height  of 
3 

0(  P)+C1  ■*  0H+C1  suggested  that  OH  vas  produced  vibrationally  hot.  This  is 
easier  to  understand  if  the  transition  state  is  non-linear.  Using  high-  .i 

level  correlated  methods,  ve  found  a  non-linear  transition  state  vith  a  C' 

CCSDT-1  barrier  of  10.2  kcal/mol  for  the  forvard  reaction  and  7.4  for  the 
reverse  [41].  Zare's  experiment  estimates  a  barrier  height  of  ~6  kcal/mol.  81 5 
-r  Jr 

K.  The  theory  and  implementation  of  analytical  gradients  for  the  full  SDTQ-  I 

MBPT(4)  [36]*  method  has  been  accomplished  this  year  for  the  first  time. 

This  vill  enable  us  to  routinely  and  efficiently  search  potential  energy  1 
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surfaces  for  minima  and  transition  states  at  the  very  attractive  and  ■■■*' 

accurate  HBPT(4)  level  of  approximation. 

L.  Similarly,  ve  have  defined  a  "response"  density  vhich  can  be  directly 
obtained  for  all  levels  of  HBPT  [36].  Using  the  response  density  ve  can 
readily  evaluate  any  one-electron  property  very  efficiently.  Ve  have  used 
this  technique  to  calculate  moments,  field  gradients  and  some  second-order 
properties  like  polarizabilities  [36]. 

M.  Using  our  recently  implemented  CCSDT  model,  the  extremely  difficult  problem 
of  the  potential  energy  curve  for  Be2  has  finally  been  resolved  [43].  Be2 

O 

has  been  found  to  have  a  minimum  near  2.5  A  vith  a  depth  of  about  2  'J 

O 

kcal/mol,  in  addition  to  the  expected  VanderVaal's  minimum  at  about  4  A  (see 
Fig.  1).  CCD,  CCSD,  and  CCSDT-1  do  not  find  this  inner  minimum  despite  the 
latter  being  within  about  1.2%  of  the  exact  result  (i.e.  full  Cl)  because 
the  1.2%  of  the  correlation  energy  accounts  for  the  small,  inner  veil.  To 
go  beyond  CCSDT-1  to  the  full  CCSDT  method,  it  is  necessary  to  fully  include 
all  effects  of  connected  triple  excitations,  T^.  In  our  recent 
implementation  of  the  full  CCSDT  model  [25],  ve  once  again  studied  the  Be2 
potential  curve.  The  full  inclusion  of  T^  accounts  for  the  inner  minimum  as 
shovn  in  the  figure. 

Another  important  consequence  of  this  study  is  that  it  offers  the  first  ab 
initio  measure  of  the  effect  of  "connected"  quadruple  excitations,  T^ 

(although  nov  our  XCC(5)  and  CCSDTQ-1  have  T^  included).  Since  Be2  is 
^studied  as  a  four  electron  system  (the  Is  electrons  are  not  correlated)  the 
difference  in  energy  betveen  the  full  Cl  result  and  CCSDT  is  the  effect  of 
T^.  This  amounts  to  about  0.3  mh  (i.e.  0.18  kcal/mol)  for  this  example. 
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N.  Continuing  our  work  using  purely  numerical  orbitals  in  CC  theory  to  attempt 
to  get  the  most  accurate  solutions  of  the  Schrddinger  equation  for  molecules 
that  have  been  obtained,  ve  studied  excited  states  of  the  dipole  bound 
anions  of  NaF,  LiCl  and  NaCl  [44].  Ve  predict  the  existence  of  such  bound 
states  at  0.012  eV,  0.009  eV  and  0.021  eV  respectively.  Ve  hope  our  study 
vill  inspire  their  experimental  observation. 

O.  Other  work  included  a  study  of  vays  to  accelerate  the  convergence  of  CC 
equations.  Results  vere  presented  for  a  series  of  small  molecules  at 
several  levels,  CCD,  CCSD,  CCSDT-1  and  the  full  CCSDT  method  [45].  Our 
reduced  linear  equation  procedure  was  shown  to  offer  the  best  convergence  of 
the  various  methods  that  have  been  used. 

P.  Ve  also  studied  the  polarizability  and  electron  affinity  of  0H~  [46]. 

Anions  have  a  very  large  polarizability  placing  extreme  demands  on  both 
basis  set  and  the  inclusion  of  electron  correlation.  To  study  the  basis  set 
effect  ve  investigated  electric  field  variant  (EFV)  orbitals,  which  are 
orbitals  that  are  explicitly  field  dependent  and  that  are  floated  off  the 
atomic  centers  to  some  optimum  location.  The  polarizability  a»47±4  a.u.  was 
almost  twice  as  large  as  the  SCF  value  attesting  to  the  importance  of 
correlation.  The  EFV  basis,  though  good,  did  not  give  much  better  results 
than  those  from  the  conventionally  chosen  basis  sets. 

Q.  Continuing  our  earlier  study  of  the  B2O  surface  in  1979,  ve  did  an  updated 
and  very  thorough  study  of  the  quartic  force  field  of  HjO  obtained  by  CCSDT- 

_1  and  other  CC/MBPT  methods  [47].  Ve  predicted  the  principal  vibrational 
frequencies  and  those  for  twenty  overtones.  The  former  vere  extremely 
accurate  2-3  cm”1  for  the  stretch  nodes,  while  the  full  set  of  frequencies 
including  the  overtones  had  an  average  error  of  less  than  about  30  cm'1.  Ve 


-8- 


also  computed  the  enharmonic  constants,  all  of  vhich  vere  vithin  2  cm' A  of 
experiment,  and  the  rotational  parameters  vhich  vere  also  exceptionally 
accurate.  Such  a  calculation  could  be  done  for  any  small  molecule  today  and 
the  results  vould  be  expected  to  be  competitive  vith  experiment  as  in  our 
H20  study. 

R.  The  IF  molecule  is  potentially  very  important  as  the  basis  for  a  laser.  Its 
dipole  moment  and  infra-red  intensity,  as  a  function  of  bond  length  is  not 
knovn.  Hence  ve  used  MBPT  to  study  the  moments  of  IF  and  other  interhalogen 
molecules  [48]. 

This  paper  also  demonstrates  the  near  satisfaction  of  the  Hellman-Feynman 
theorem  for  high-order  MBPT  vave functions,  vhich  is  an  important  result 
contributing  tovard  our  development  of  the  nev  class  of  UCC  methods 
discussed  above. 

S.  Additional  applications  of  our  CC/MBPT  method  focused  on  the  activation 
barrier  and  transition  state  for  the  autoisomerization  of  cyclobutadiene 
[16]  and  the  isomerization  energies  and  related  properties  of  the 
cyanoborane  complexes,  BH^CN-,  BH^NC',  HCNBHV  and  HNCBH^  [49]. 

T.  Formulae  to  evaluate  the  spin  multiplicity  in  unrestricted  Hartree-Fock 
(UHF)  based  CC/MBPT  methods  vere  derived  and  applied  to  demonstrate  hov  CC 
theory  eliminates  most  of  the  spin  contamination  in  such  a  calculation, 
enabling  UHF/CC  theory  to  provide  good  potential  energy  surfaces  near  atomic 

^dissociation  limits  [50]. 
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